Phase transition of the one-dimensional coagulation-production process 
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Recently an exact solution has been found Q for the Id coagulation production process: 2A — > A, 
A0A — > 3A with equal diffusion and coagulation rates. This model evolves into the inactive phase 
independently of the production rate with t" 1 / 2 density decay law. Here I show that cluster mean- 
field approximations and Monte Carlo simulations predict a continuous phase transition for higher 
diffusion/coagulation rates as considered in [jjj. Numerical evidence is given that the phase transition 
universality agrees with that of the annihilation-fission model with low diffusions. 
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One-dimensional, non-equilibrium phase transitions 
have been found to belong to a few universality classes, 
the most robust of them is the directed percolation (DP) 
class [|]||. According to the hypothesis of all con- 
tinuous phase transitions to a single absorbing state in 
homogeneous systems with short ranged interactions be- 
long to this class provided there is no additional symme- 
try and quenched randomness present. 

Recent studies on the annihilation fission (AF or 
PCPD) process 2 A -> 0, 2 A -> 3 A §-§ found evidence 
that there is a phase transition in this model that does 
not belong to any known universality classes. This model 
without the diffusion of single particles was introduced 
originally by . The renormalization group analysis of 
the corresponding bosonic field theory was given by ||. 
This study predicted a non-DP class transition, but it 
could not tell to which universality class this transition 
really belongs. An explanation based on symmetry ar- 
guments are still missing but numerical simulations sug- 
gest p|)p[ that the behavior of this system can be well 
described (at least for strong diffusion) by coupled sub- 
systems: single particles performing annihilating random 
walk coupled to pairs (B) following DP process: B — > 2B, 
B — > 0. The system has two non-symmetric absorbing 
states: one is completely empty, in the other a single 
particle walks randomly. Owing to this fluctuating ab- 
sorbing state this model does not oppose the conditions 
of the DP hypothesis. Some exponents are close to those 
of the PC class (^|| but the order parameter exponent 
(/3) has been found to be very far away from both of the 
DP and PC class values Q]. In fact this system does not 
exhibit neither a Zi symmetry nor a parity conservation 
that appear in models with PC class transition. It is 
conjectured that this kind of phase transition appears 
in models where (i) solitary particles diffuse, (ii) particle 
creation requires two particles and (iii) particle removal 
requires at least two particles to meet. 

In this paper the following one-dimensional coagula- 
tion production processes will be investigated: 

a) Spatially symmetric coagulation production pro- 
cesses: 



A0A 



3A, 



2A ^ A0, 2A ^ 0A, 
A0 0A 



(1) 
(2) 
(3) 



b) Spatially asymmetric coagulation production pro- 
cesses: 



AA0 m> 3A, 0AA 3A, 
2A ^ A0, 2A ^ 0A, 
A0 A 0A 



(4) 
(5) 
(6) 



Both versions fulfil conditions (x— iii . ) but Henkel et al. 
|jj show that for d = c the symmetric version always 
evolve into the inactive state with p oc t~ 0,5 scaling 
law. They argue that the asymmetric version displays a 
non-equilibrium phase transition. The difference is said 
to be similar to the hard-core effects observed in one- 
dimensional models |T^-|T6[]. Hard-core particle exclu- 
sion effects can really change both the dynamic (lj 
and static Jl5|,|l6| behavior of one dimensional systems by 
introducing blockades into the particle dynamics but in 
this work I argue that not this kind of hard-core effects 
responsible for the lack of phase transition. 

One can quickly check by simulations that for d < c 
the density in the asymmetric version decays in much 
the same way - with p cx t~ - 5 scaling law - as in case 
of the symmetric version. Furthermore I shall show that 
if the coagulation rate is smaller than the diffusion rate 
particles can escape before removal an active phase will 
emerge with a continuous phase transition belonging to 
same class that was found in the AF model for weak dif- 
fusion. Therefore both version exhibit qualitatively the 
same phase diagram. 

To prove this first I shall apply cluster mean-field ap- 
proximations (GMF) Jl7],[l8| , which can predict phase di- 
agrams qualitatively well. The mean-field equation for 
the steady state of both version is 



= f(l- PA ) P 2 A 



CPA 



(7) 



1 



where pa is the probability of As at a given site. Note 
that the diffusion rate d does not play a role in this ap- 
proximation. By introducing the parametrization c = 
p(l — d), f = (1 — p)(l — d) - that is similar to that of 
the PCPD model - this has the solution: 



PA 



2p-l 

v-i 



for p < 1/2 and p = if p > 1/2. Therefore an active 
state appears in the mean-field approximation already. 

For higher order cluster mean-field approximations 
similar scenario can be found, but one has to treat the 
two versions separately. The density in pair approxima- 
tion for the symmetric version is: 



(p - 1) p 2 - 2 dp (p 2 + 2p - 2) 



d 2 (p 3 



bp 2 + Ap - 4) 



0-1) p 2 - 2 dp (p 2 + 2p- 2) + d 2 (p 3 + 5p 2 - 4) 

(9) 

One can easily prove that if the coagulation rate is equal 
to the diffusion rate d = p(l — d)/2 this gives a single 
p = absorbing state solution in agreement with |IJ] . The 
steady state solution with positive density is possible if 



d> 



P -P 
' + 3p - 2 



(10) 



This gives the phase boundary in pair approximation that 
is a continuous unlike in case of the PCPD model § (see 
Fig.0). 
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FIG. 1. Phase diagram of the symmetric coagula- 
tion-production model. Dotted line : mean-field approxima- 
tion, dashed line : pair approximation, squares: simulation 
results. The circles show d c as the function of c. Lines con- 
necting symbols are used to guide eyes only. 

The pair density in this approximation 

((p-i)p-d (p 2 + 3 P -2)) 2 (p-iy 1 



(p 



1) p 2 -2dp (p 2 + 2p- 2) + d 2 (p 3 + 5p 2 -4) 

(11) 



has a leading order singularity all along the phase tran- 
sition line 



c oc (p c - pf 



(12) 



suggesting one universality class unlike in the case of 
PCPD model 0. 

The GMF solutions for N = 3, 4, 5, 6, 7 block sizes have 
been determined numerically at d = 0.2. The approxi- 
mation level is constarined by the numerical stability of 
the fixed point solution in the multi-dimensional space 
of iV-block probability variables. As Figure § shows the 
Pn density curves of different approximations converge 
to the simulation results. 
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FIG. 2. Cluster mean-field approximations in the symmet- 
ric coagulation-production model for d = 0.2. The curves 
correspond to steady state density solutions as the function 
of p, for iV = 2, 3, 4, 5, 6, 7 (right to left). The circles with er- 
ror bars represent simulation results. Insect: Corresponding 
coherent anomaly amplitudes with a power-law fitting. 

Using these data an estimate can be given for the or- 
der parameter density exponent p oc \p — p c \$ using the 
Coherent anomaly method (CAM) [^9), which has been 
proven to give precise estimates for the DP |2(J and PC 
PH classes. According to CAM the amplitudes a(N) of 
the cluster mean-field singularities scale in such a way 
that 



a(N)oc\p c (N)-p c f-P* 



(13) 



the exponent of true singular behavior can be estimated. 
From the mean-field solution (^) one read-off that 0m F = 
1. The critical point p c can be estimated either by ex- 
trapolating on the GMF results or by simulations. Lin- 
ear extrapolation at d = 0.2 for p c (l/N — > 0) gives: 
p c = 0.182(2). Monte Carlo simulations on large sys- 
tems - discussed below - give a more precise estimate: 
p c = 0.17975(8). The amplitudes a(N) near p c (N) are 
determined by linear fitting from the pn{p) data and 
shown in insect of Fig.||. as the function of p c (N) —p c - A 
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power-law with exponent ft — (3mf = —0.43(3) can fairly 
well applied for points corresponding to N > 2 approx- 
imations giving an estimate: ft = 0.57(3), which agrees 
well with former results for the AF model with small dif- 
fusion rates ||. 

Monte-Carlo simulations of the symmetric process 
started from fully occupied lattices of size L — 40000 
show a phase transition for d = 0.2 and p c = 0.17975(10) 
(see Fig.||). The local slopes of the density decay: 



eff(t) 



\n[p(t)/p(t/m)} 
ln(m) 



(14) 




FIG. 4. Effective order parameter exponent results. Linear 
extrapolation results in j3 = 0.57(1). 



(where we use m = 8 usually) at the critical point go to 
exponent a by a straight line, while in sub(super)-critical 
cases they veer down(up) respectively. 
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FIG. 3. Local slopes of the density decay in the symmetric 
coagulation production process. Different curves correspond 
to p = 0.1795, 0.1797, 0.1798, 0.1799 0.18 (from bottom to 
top). Throughout the whole paper t is measured in units of 
Monte-Carlo sweeps (MCS). 

For the critical point (p c — 0.17975(8)) one can esti- 
mate that the effective exponent tends to a = 0.263(9), 
which agrees with results for the AF model |^|j9| again. 
For other d-s similar results have been found. 

In the supercritical region the steady states have been 
determined for different e = p— p c values. Following level- 
off the densities were averaged over 10 4 MCS and 1000 
samples. By looking at the effective exponent defined as 



Peff^i) = — ; ; , (15) 

lnei - lne 4 _i 



one can read-off: (3 e ff — ► ft — 0.57(1), which is in good 
agreement with the exponent of the AF model for weak 
diffusion determined by coherent anomaly method and 
simulations @. 

The simulations and the cluster mean-field approxima- 
tions show that if the diffusion rate is lowered this phase 
transition disappears and the system will decay with the 
p cx t~ - 5 law independently of / in both versions. As 
expected the asymmetric version exhibits a phase transi- 
tion with the same universal properties as the symmetric 
version. Example for d = 0.2 the transition point is at 
c = 0.359(1), f/2 = 0.4409(1) with the decay exponent 
a = 0.27(1). 

In conclusion coagulation production models exhibit a 
phase transition if the diffusion is fast enough. The spa- 
tial symmetry of the production process has been found 
to be irrelevant as in case of the AF process |1 . The crit- 
ical behavior agrees well with that of the AF model in its 
weak diffusion rate region. An open question is that why 
can one not see the cyclically coupled behavior in this 
model similarly as in the PCPD model as d — » 1. The 
corrections to scaling are getting very strong in this limit 
that make numerical solutions very confusing, but one 
has to realize that the B — > process of pairs (present 
in AF) is missing in this model. Therefore a single uni- 
versality class in this model and two distinct classes in 
the AF model is likely. This conjecture is strengthened 
by the pair mean-field results : one obtains analytically 
the same singular behavior here and two distinct singu- 
lar behavior in case of the PCPD model along the phase 
transition line. 



TABLE I. Summary of results 



d 


0.1 


0.2 


0.4 


0.7 


Pc 


0.1129(1) 


0.17975(8) 


0.2647(1) 


0.3528(2) 


a 




0.263(9) 


0.268(8) 


0.275(8) 


13 




0.57(1) 


0.58(1) 


0.57(1) 


/3cAM 




0.57(3) 
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